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SUMMARY 


The free turbulent mixing of a supersonic Jet of Mach number 1.6 has 
been experimentally investigated. An interferometer, of which a descrip- 
tion is given, was used for the investigation. Density and velocity 
distributions through the mixing zone have been obtained. It was found 
that there was s imil arity in distribution at the cross sections investi- 
gated and that , in the subsonic portion of the mixing zone , the velocity 
distribution fitted the theoretical distribution for incompressible flow. 
It was found that the rates of spread of the mixing zone both into the jet 
and into the ambient air were less than those of subsonic Jets. 


INTRODUCTION 


A considerable amount of work, both theoretical and experimental, 
has been done in the past on free jets. Most of the theoretical work 
has been based on Prandtl's concept of the "mixing length." Tollmien 
(reference l) and Gortler (reference 2) treated the turbulent mixing of 
incompressible jets. Abramovich (reference 3) published a theory of 
the free subsonic jet in which effects of compressibility were included. 
Experimental results on velocity distribution in the mixing region and 
rate of spread of the mixing region both into the jet and into the 
surrounding air have been published by various investigators (references If 
to 8). Most of the experimental work has dealt with constant-density 
cases j that is, the density of the jet was the same as that of the 
ambient air with which it mixed, and the velocity of the jet was quite 
small. The present paper reports results of an experimental investiga- 
tion of the mixing of a supersonic jet. In this jet the Mach number was 
1.6 and the density was about one and one half times the density of the 
ambient air. Measurements were made with an interferometer of the density 
variation across the mixing zone in the region near the nozzle. From the 
density variations the velocity variations have been calculated. A full 
description of the interferometer is given, together with discussions of 
the technique of adjusting the interferometer, the theory of fringe 
formation, and the method of evaluating interferograms . 
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SYMBOLS 


a 

"b 

c 


constant 
fringe spacing; constant 
constant 
constant in Snell's law 


icing, cou 


c 2 


C = 

k 

l. 

L 

n 

XI 1 


constant of integration 

specific heat of air at constant pressure 

L(n - 1) 

Gladstone -Dale constant 
mixing length 

distance throu^i test section 

index of refraction of air of density p 

index of refraction of air of density p' 



S(y) 


T 

m 

stag 

u 

^o 

V 


dimensionless fringe shift, Y/b , where density is p' 
static temperature of air 
stagnation temperature of air 

Telocity of air 

velocity of air in free stream. 

velocity of light 
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V* 


x 

2 
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Y 

z 

X 

X o 

X’ 


velocity of light in vacuum 

velocity of light in air of density p* 
coordinate 

retardation of ligjit "beam that causes fringe shift Y 

coordinate 

fringe shift 

coordinate 

wavelength of light 

wavelength of light in vacuum 

wavelength of light in air of density p' 


I = p + qy 


P>P' > p atm density of air 


or 


4 


scale factor determined hy comparing experimental 
theoretical velocity distributions 

angle of incidence of light ray 


APPARATUS 

Jet 


The open jet is operated hy the air from a 500-cubic-foot storage t ank 
in which the initial pressure is about l 80 pounds per square inch and from 
which a pipe leads to a supersonic nozzle. The exit end of the nozzle is 
open to the atmosphere. In the pipe that runs from the tank to the nozzle 
is a hand-operated gate valve , an air-operated quick-opening valve , and 
a pressure regulator. The purpose of the pressure regulator is to reduce 
the pressure of the air going into the nozzle to a value that will result 
in the pressure of the air } as it leaves the nozzle } being atmospheric. 

Th© result is that no strong shock waves or expansion waves originate 
at the rim of the nozzle to adjust the jet pressure to atmospheric . 

The results reported herein were obtained with a nozzle that gave 
a Mach number of 1.6. The nozzle was of the two-dimensional type and 
was constructed of steel. It had an exit opening of 3 by 3 inches. The 
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coordinates of the nozzle "blocks were calculated "by the method of charac- 
teristics. No allowance was made for the thickness of- the boundary layer, 
and the nozzle gave the Mach number for which it was designed. A photo- 
graph of similar nozzles is shown as figure 1, together with the section 
for the transition from a 6-inch circular pipe to the rectangular entrance 
to the nozzle. 


Interferometer 

Introduction .- The Mach-Zehnder type of interferometer was originated 
by L. Mach and Zehnder for use in studying phenomena of gas dynamics. 

The instrument iB useful in the study of gas-flow problems because it 
gives an instantaneous record freon which can be calculated the variation 
of gas density throughout a flow field. It is particularly applicable to 
the study of high-speed gas flows, because, by using light waves, it makes 
unnecessary the insertion into the flow of probes or other measuring 
instruments that would disturb the flow. The interferometer gives quanti- 
tative results in which a rather high degree of accuracy can be obtained. 

It has a sufficient range of sensitivity to measure both small density 
changes, as in a weak Prandtl-Meyer expansion, and large density changes, 
as across a strong shock wave. 

The interferometer was used, but not intensively, by Mach and Zehnder. 
It was applied to the study of subsonic aerodynamics by Zobel (reference 9) 
Since then it has been applied by Ladenburg and his coworkers to the study 
of phenomena in supersonic flow (reference 10 )• 

J 

The basic arran ge ment of the Mach-Zehnder interferometer is shown in 
figure 2 . Light from a source S is made into a beam of parallel rays 
by a collimating lens system L]_. This beam falls on the splitter plate 

S-^ where it is split into two beams. Part of the original beam of light 
is reflected by the splitter plate S-|_ and part is transmitted by it. 

The part that is transmitted goes to the mirror M-j_ where it is reflected 
onto the splitter plate Sg. A portion of the beam is transmitted by Sg 
and is not used. The other portion is reflected by S2, passes through 
the lens Lg, and f all s on a screen or a photographic plate P. 

The light that was reflected by Sq_ likewise goes to a totally 
reflecting mirror Mg from which it is reflected onto the splitter plate 
Sg. At S 2 a part of the beam is reflected and is not used, and the 
remainder is transmitted, passes through the lens Lg, and falls on the 

screen or the photographic plate P. This arrangement fulfills one of 
the necessary conditions for -the interference of the waves in two beams 
of light y that is, that the two beams originate in the Bame light source. 
From a practical standpoint the arrangement also permits the condition to 
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"be met that the two beams be widely enough separated in space that the 
disturbance to be studied can be introduced into one of the beams without 
disturb ing the other. The disturbance , or the "test section," can, of 
course, be located anywhere in either of the two beams. In the apparatus 
described in the present paper the test section was situated midway 
between the mirror ^ and the splitter plate Sg » 

Theory of ideal fringe formation .- The two beams of light that reach 
the photographic plate P appear to come from separate sources that are 
situated somewhere to the left of the mirror H,. By proper orientation 

of the two splitter plates and the two mirrors, the two beams of light can 
be made to appear to cro3S each other, as is shown in figure 3* If, for 
the moment, it is assumed that each beam is composed of strictly parallel 
and monochromatic rays , then the ware fronts can be represented by equally 
spaced straight lines perpendicular to the direction of propagation, as is 
also shown in figure 3* Here the straight lines represent the "crests" of 
the waves. Midway between two successive "crests" are the "troughs" of the 
waves. Where two lines intersect, two crests occupy the same position in 
space, reinforce each other, and cause an increase in the amplitude of 
vibration and an increase in the intensity of the light. Where a line 
intersects a point that is midway between two adjacent lines of the other 
beam, a crest and a trough interfere destructively, so that a decrease in 
the amplitude of motion and a decrease in the intensity of the light results. 
The plate P will therefore be crossed by parallel, horizontal (in this case 
lines of alternately weak and strong intensity. These lines are the 
interference fringes. They can be oriented in any direction by proper 
rotation of the two splitter plates and the two mirrors about two axes, 
one in a plane parallel to that of the paper and one perpendicular to the 
plane of the paper. 

If now a disturbance is produced in the air that is in the test section, 
and this disturbance changes the speed of the light in the beam that actually 
traverses the test section, then the wave fronts in that beam will be 
advanced or retarded, and the positions of the fringes will be shifted down 
or up. Prom the shift in the fringes the average speed of the light and, 
from that, the average density of the air can be calculated. The only gas- 
flow caseB for which the density field can be evaluated are the one- 
dimensional case of uniform density throughout the field, the two-dimensional 
case of uniform densities along- lines (or planes) that lie parallel to the 
light-beam direction of propagation, and the three-dimensional case of 
axially symmetric densities, such as the flow field about bodies of 
revolution at zero angle of attack. 

Theory of -practical fringe formation .- So far it has been assumed 
that th9 light is strictly monochromatic and is a beam of parallel rays. 

Bat parallel light can be obtained only from a truly point source. 

Because in practice neither monochromatic light nor a point source is 
used (because, with available means of approximating these to a very high 
degree, the intensity of the light would be too .small), neither condition 
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is actually met. As has been pointed out before (reference 10 ), the 
complete theory of fringe formation for the actual case has not yet been 
given. Schardin, however, has shown (reference ll) that the following is 
the case. With a point source and a strictly parallel beam, fringes are 
formed at ai 1 points along the two axes I]_ - I]_' and I2 - 12* , where 

the two beams overlap. (See fig. 4.) Bat in an actual case, in which a 
light source of finite extent is used, the lines 1^ - I-^' and I2 - 

become very numerous. The most distinct fringes are then formed where the 
various lines between 1^ 1^' and between I2 and Iq' intersect. 

The interferometer should be so adjusted, by rotation of the splitter 
plates, that this intersection is in the center of the test section. 

For obtaining the greatest contrast between fringes , it is also 
necessary that the optical-path length through the interferometer be the 
same for the two beams. In other words, it is necessary that the optical 
path for the two beams be nearly the same from the time the original beam 
is split into two, at the splitter plate until the two are reunited 

at the splitter plate This condition becomes the more Important the 

more the light departs from being monochromatic . For a given setting of 
the splitter plates and the mirrors, the fringes have a given spacing for 
each wavelength. The greater the wavelength, the greater the spacing 
between wavefronts and consequently between fringes. Thus, fringes 
produced by a large numb er of wavelengths coincide only at one point , the 
center of the band of fringes^ On either side of this center the fringes 
get out of step with each other and the contrast between light and dark 
becomes less ainrl less. With the optical paths equal, the center of the 
band of fringes, where the contrast is greatest, lies at the center of 
the test section. 

Description of interferometer .- The interferometer with which the 
results reported in the present paper were obtained was designed and 
constructed at the Langley Laboratory and is installed in the boundary 
layer laboratory of the Physical Research Division. The base on which 
the splitter plates and mirrors are mounted is a one-piece iron casting 
in the foim of a four-leaf clover. (See figs. 5 ^d 6.) The assembly 
is supported in a vertical plane at its center by a single mount , which 
is attached to a framework of structural steel that is welded to a steel 
table. The table rests on steel plates that are bolted to the concrete 
floor of the second story of the building. The building happens to house 
numerous motors, compressors, and other sources of vibration. It is within 
50 feet of a projectile gallery, a 500 -horsepower wind tunnel, and a 
1000 -horsepower wind tunnel. ALL of these pieces of equipment cause 
vibrations of considerable amplitude in the second-story floor, which is 
supported by columns only every 20 feet or so. 
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Interferometers are quite delicate instruments, as far as adjustment 
is concerned, and are quite subject to being thrown out of adjustment by 
mechanical vibrations and by temperature changes. When the present 
instrument was designed it was hoped to produce an instrument that would 
not cause an inordinate waste of time in keeping it in adjustment. The 
result has been quite satisfactory. Apparently because of its symmetrical 
design (and probably also because cast iron has a low coefficient of 
thermal expansion) , no effects of temperature changes on the adjustment of 
the interferometer have ever been noticed. Ordinarily the interferometer 
is enclosed in a sheet-aluminum case. Even when it is used with one side 
of the case removed and the temperature in the room is changed by opening 
the windows, or when the operators stand close to the interferometer, no 
change in the fringes is noticed. 

In an effort to keep the interferometer in adjustment despite vibrations 
it was at first suspended from a framework by springs that gave the system 
a natural frequency of about 1 cycle per second. It was soon found, when 
studying the boundary layer on a flat plate, that the beam of li^it could 
not be kept lined up parallel' with the flat plate. Moreover, at that time 
the glass plates and mirrors were supported at four points, and they could 
not be clamped very tightly without causing strains in the glass and 
distortion of the fringes. As a result, the interferometer would not stay 
long in adjustment, despite the fact that it was suspended by springs. 
Accordingly, the supports of the plates and mirrors were changed to the 
three-point type (see fig. 7) and were clamped down very tightly without 
distorting the fringes. The spring support was discarded and the inter- 
ferometer table was placed on the floor as just described. As a result, 
the interferometer has stayed In adjustment for some length of time. 

In fact, in the past 9 months no adjustments of any kind have been made, 
and the interferometer has stayed in perfect adjustment. 

Each of the two plates and two mirrors is rotatable about two axes, 
one in a plane parallel to that of the interferometer base (fig. 6) and the 
other perpendicular to the plane of the interferometer base. For. rotation 
about the first of these axes, the mirror holder is mounted in Journal 
bearings. (See fig. 7.) These bearings are clamped tight, Just short of 
binding. A lever is attached to the axle. The end3 of micrometers press 
against the lever. Rotation of the plate or mirror is effected by 
advancing one micrometer while retracting the other. When the proper 
position is reached, both micrometers are screwed very tightly against 
the lever. For rotation about the other axis, the plate or mir ror housing 
rotates on ball bearings. Again a lever is moved by means of two micrometers. 
(The micrometers provide an economical method of obtaining precision screw 
threads. The scales on the micrometers are not used.) 

For adjustment of the white-light fringes; that is, for adjusting 
the difference in optical-path length of the two beams, where optical- 
path length is defined as the integral of the product of index of 
refraction and the differential of geometric path length, an adjustment 
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was provided whereby the entire housing of mirror M-j_ could he translated 

in a direction perpendicular to the plane of the mirror. The ways on 
which the hous ing moved were found, however, to he not sufficiently smooth, 
and when the mirror was moved, the plane of the mirror changed, and thus 
the orientation of the fringes changed. The screw was also found to he too 
coarse. Therefore, for easier adjustment of the white-light fringes, a 
compensating plate was Installed in each "beam. (C^ and C2 , fig. 6) • 

Each plate is rotatahle about two. mutually perpendicular horizontal axes . 

One plate is placed in a horizontal position, and the white-light fringes 
are adjusted hy rotating the other plate. The fact that the path in glass 
may he different for the two beams does not matter, because the light is 
too nearly monochromatic for dispersion to have an appreciable effect. 

The splitter plates and the mirrors are 4 inches square and l/2 inch 
thick. The plates are polished flat on each side to a tolerance of l/lO 
wavelength. The angle of the wedge foxmed by the two sides of a plate does 
not exceed 2 seconds of arc. One side of each plate is coated with zinc 
sulphide for increased reflection. Use of this coating is a departure 
from the conventional method. This coating has the advantage over the 
conventional one, which is a thin coating of silver or aluminum, in that 
no light is absorbed by the coating. The ideal coating would transmit 
50 percent of the light and reflect 50 percent. Since each beam is 
transmitted once and reflected once at splitter plates, the net result would 
be two beams going into the camera, each of which had 25 percent of the 
intensity of the original beam. For the zinc sulphide coating that was 
used, approximately 40 percent is reflected and 60 percent transmitted. 

This is satisfactory, however, because each beam into the camera then has 
24 percent of the original intensity (neglecting losses at the mirrors). 

The two mirrors are flat to a tolerance of l/lO wavelength and are 
front- surfaced with rhodium. Because of the low reflectance of rhodium 
(about 75 percent), it is expected that aluminum would be more satisfactory. 

Description of light source and camera .- A satisfactory light source 
for application of interferometry to the study of flow phenomena must 
provide nearly monochromatic light in a nearly parallel beam.. An 
additional requirement for taking inter? ero grams of hi^i-speed flow is 
that the duration of the light be sufficiently short. For taking 
interfere grams that are free from blurring of phenomena in an open, 
or free, supersonic jet, exposure times of the order of three microseconds 
or less must be used, (inherent in free supersonic Jets are high-frequency 
vibrations.) Because meeting each of the requirements of being mono- 
chromatic, parallel, and of short duration tends to reduce the intensity, 
one must add the fourth and obvious requirement of sufficient intensity. 

The light-source problem has been solved by the same general method 
that was used by the Princeton group (reference 10 ), that is, by use of a 
high-voltage magnesium spark and a monochromator. A diagram of the 
complete light-source optical system is shown in figure 8. The two 
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electrodes are magnesium rods of i-inch diameter with rounded tips 

approximately 1 inch apart. Each magnesium rod is held in a concentric 
hollow "brass rod. A press fit is used for mechanical strength, and where 
the magnesium rod enters the brass rod, a weld is made. Each "brass rod 
is soldered to a thin copper sheet, which extends for about an inch and 
is then soldered to a terminal on each of two condensers. Each of the 
two condensers, which are thus connected in parallel, is a 15, 000-volt 
pyranol condenser and has, according to the manufacturer’s specifications, 
a capacitance of 1.0 microfarad and an inductance of 0.5 microhenry. The 
condensers are charged, through a 10 megohm isolating resistor, to 
16,000 volts. The discharge is initiated by means of a tickler spark 
between one of the magnesium electrodes and an auxiliary magnesium 
electrode placed between the other two and slightly to one side. The 
tickler spark is produced by the potential across the secondary of an 
automobile ignition coil that occurs when the 6-volt circuit through the 
primary is broken. This arrangement gives satisfactory control over the 
main discharge and facilitates synchronizing it with the opening of the 
camera shutter . The duration of the light from the main discharge was 
measured with a rotating-mirror apparatus and was found to be of the 
order of 3 microseconds. 


The monochromator was constructed from a Bausch and Lomb "Simplified 
Constant Deviation Prism Type" Laboratory Spectrometer. The collimator 
objective, of aperature f/8, and the entrance slit of the original 
spectrometer were removed; and were replaced by a system of greater 
aperture that was constructed from available lenses. The lens 

(fig. 8) is an f/l.6 Kodak Anastigmat of 50 -millimeter focal length. 

This lens focuses the light from the spark discharge on the adjustable 
entrance slit. (No pinhole is used at the spark.) The light is then 
made into a parallel beam by the lens L4, which is an f/2.3 Bausch and 

Lomb B altar of 2 -inch focal length. The parallel beam then goes through 
the constant-deviation prism of the Pell in-Broca type. The beam is then 
focused on the adjustable exit slit by the lens L5, which is an f/8 lens 
of 6-inch focal length and is the original lens furnished with the 
spectrometer. The exit slit is at the focus of the mirror Lj_. The 

light from the slit is turned at right angles by a small prism, which 
is placed slightly off the axis of L^. The mirror L-l is a parabolic 

mirror of 4-inch diameter and 36-inch focal length and is front- surfaced 
with chrome - aluminum • 

It is realized that considerably more light could be obtained from 
the system if the f-numbers of the various lenses and the mirror were 
properly matched, in order that the image of the entrance slit, which 
falls on the exit slit, would be more nearly the same size as the 
entrance and the exit slits# The system as used at present, however, works 
satisfactorily , and therefore no changes in it are contemplated* 
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The light from the. green triplet of magnesium is used, which has 
wavelengths of 5167 , 5173; and- 51 8 ^ angstrom units. Of course, at the 
high exciting voltage that is used, the lig£it does not consist only of 
these three wavelengths "but is nearly a continuous spectrum. It was 
found that the light from this region produced much more satisfactory 
interference fringes than that from any other region of the spectrum of 
a high -voltage magnesium spark. (The Princeton group reported (reference 10 ) 
that they obtained best results by using the blue line at 4h8l angstroms.) 

The monochromator was set at 5170 angstroms, and the slits were set at a 
width of 0.3 millimeter. A band about 30 angstroms wide was passed by the 
exit slit, and about 180 usable fringes were obtained. It was found 
that the appearance, or contrast, of the fringes could be improved by 
reducing the length of the exit slit to 3/32 inch. 

The fringes were photographed with an Eastman Anastigmat aerial 
camera lens of 13* 5" inch focal length and f/3*5 aperture. Kodak Lina- 
graph Ortho film was used. The Linagraph Ortho is a very fast ortho- 
chromatic film of moderate contrast and hi^i resolving power and is 
designed for photographing high-speed transient phenomena on green- 

fluorescing cathode-ray screens. Negatives of about 1^-inch diameter 

were taken. (Inasmuch as the film was used in the 35-mm size, a portion 
of the light did not hit the film. Since the interf erograms shown in 
the present paper were taken, the film has been changed to the 70 -mm 
size.) The moderate grain size of the f ilm permitted enlargements of 
sixteen or more times the diameter of the negative, or seven or more 
times actual size. 

A Kodak Supermatic shutter was placed at the focus of the camera 
lens. The shutter was set at l/50 second and was synchronized with the 
light-source spark. For convenience in moving the film, a slightly 
modified Argus C-3 camera, with the lens removed, was used as a film 
holder. 

The procedure that was followed for each interf erogram of the flow 
was to take first an interf erogram with no flow, then one with a 
transparent ruler in the test section for establishing the scale , then 
one with flow, then a final one without flow. The film was developed 
in Kodak D -72 for the exceptionally long time of lh minutes at 68 ° F. 

Adjustment of the interferometer .- When the interferometer is first 
sob up, a number of adjustments must be made in order to obtain fringes 
and to orient them properly. Most of the adjustment procedure that was 
used is a more -or- less standard procedure , but part is original. 

The first step in the initial adjustment of the interferometer is 
to make the reflecting surfaces of the two splitter plates and the two 
mirrors nearly parallel. This is done by making the two mirrors as 
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nearly parallel as possible by eye, and then leaving them there, because 
it is possible to produce all fringe orientations by adjustment of the 
splitter plates alone (except for adjustment for white-light fringes which 
is made with the compensating plates). Then the plates are set nearly 
parallel to the mirrors, by simple inspection. A small light source is 
placed 15 feet or more away from the interferometer, and the light is 
directed at plate S-^. Two sets of cross hairs are set up, one near the 
light source and the other near plate Sj_. Two screens are set up on 
the opposite side of the interferometer in such positions that a lens 
placed after plate Sg focuses one set of cross hairs on one screen and 
the other set on the other screen. Each set of cross hairs produces two 
images on one screen.. One image is produced by the lig£it that is trans- 
mitted throu gh the lower part of the interferometer, and the other image 
by the light that goes through the upper part. The splitter plates are 
then so adjusted that the two images of one set of cross hairs coincide, 
or merge into a single image. This procedure is repeated for the other 
set of cross hairs. As the adjustments for the two sets of images are 
not In dependent, they must be continued until the images of both sets of 
cross hairs appear to be single. The light source is now replaced by a 
monochromatic light source. This light source should be sufficiently 
close to being monochromatic that many hundreds of fringes can be 
produced. The reason for this is that at this stage of the adjustment 
the two optical paths through the interferometer may be considerably 
different. If such is the case, and a light source that would give 
relatively few fringes is used, then the fringes cannot be found 
because their apparent location would be above or below the test section 
and out of the light beam. A satisfactory light source is a sodium-arc 
lamp, or a General Electric B-H6 mercury lamp operated at an undervoltage 
of about 60 volts. The lens in the beam emergent from the interferometer is 
replaced by a telescope placed some 15 feet beyond the interferometer. 
Interference fringes generally can then be brought into focus at some point 
between the light source and the telescope. If they cannot be found, either 
the elimination of double cross-hair images has not been sufficiently achieved 
or the difference in optical -path length is too great. Inspection of each 
cross-hair image (on the screens) with a magnifying glass will usually reveal 
that the elimination of doubling has not been completely effected. 

Once the fringes have been located, the monochromator is set up. 

A B-H6 mercury lamp is placed at the entrance to the monochromator. 

The slits are narrowed to about 0.3 millimeter and the prison set to 
pass the 5^-60-angstrom line. The next adjustment is to tilt the fringes 
into the desired orientation, say a horizontal position, by rotating a 
splitter plate. 

The next two adjustments are to move the fringes into the test section 
and to adjust them to the desired width, or spacing. These two adjustments 
are made simultaneously by alternate rotation of the two splitter plates 
and by observing through the telescope the. location and the spacing of 
the fringes. This adjustment is not difficult . For example , suppose 
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the fringes are located "between S2 and the telescope and are too narrow* 

The first step is to determine which of the two beams goes through which 
path in the interferometer. This is determined by focusing the telescope 
on the light source and blocking off one of the paths. It is desired to 
have the beams cross in the test section , midway between Mg and £>2. 

If the slope of the beam that passes through the test section is positive 
with respect to the slope of the other beam, as shown in figure 4, then 
plate Sg is rotated counterclockwise to bring the fringes back almost 
to the test section. This adjustment may make the fringes too broad. Then 
S-l is rotated clockwise to move the fringes the rest of the way to the 

test section and at the same time to narrow them. 

The next adjustment is to make the two optical-path lengths through 
the interferometer equal. With the fringes in focus in or near the test 
section, the slits are opened until the fringes visible in the telescope 
become faint. By rotating a compensating plate, the fringes are caused to 
appear to pass vertically through the test section until they disappear. 

The position of the compensating plate ,is noted, and the plate is rotated 
in the opposite direction until the fringes with the greatest contrast 
have passed through the test section and the fringes again disappear. The 
compensating plate is then positioned halfway between the two positions where 
the fringes disappeared. The interferometer is then nearly in adjustment 
for white-light fringes. Either the monochromator is removed and a source 
of white light is used, or the white light is placed Just ahead of the 
first splitter plate without disturbing the monochromator. (A convenient 
point source is the zirconium-arc light. ) Then a slight adjustment of a 
compensating plate will bring the white-light fringes into view. If the 
disturbance of which an interferogram is to be taken cont ains a region in 
.which there is a density gradient of considerable magnitude, such as a 
boundary layer, the fringes will be crowded together in that region, and 
distinguishing individual fringes may be difficult. It is advisable to 
place the white-light fringes in such a position that they will move, 
when the disturbance is produced, into the region of density gradient 
and thereby provide the fringes of greatest contrast in that region. 

The final adjustment is to remove what might be called "twist" 
from the two beams. For horizontal fringes the two splitter plates have 
been rotated about horizontal axes and the two images of the light source, 
as seen in the telescope, lie one above the other. It may be, however, 
that the two beams do not lie in the same vertical plane. If they do not, 
sharp fringes can be observed only when the source is a line of zero 
width (and is also vertical, for horizontal fringes). As it is necessary 
that the light source have a finite width, in order to give enough light, 
then the fringes will appear considerably blurred unless the rotation of 
the plates is corrected in- order that the two beams lie in the same 
vertical plane. To check this, the telescope is removed and the eye 
placed in the emergent beam some distance away from the interferometer. 
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The fringes may then appear to be cocked at an angle to the horizontal. 

If so, then as one walks toward the interferometer the fringes will 
slowly rotate back to the horizontal position which they had when viewed 
with the telescope focused on the test Section. By alternately rotating 
the two plates about their non-horizontal axes, the two beams can be 
swung until they both lie in the same vertical plane and the fringes 
appear horizontal when viewed from both a close and a far position. 

When all these adjustments have been made , the following will be true : 
. (a) The fringes are centered in the test section. 

(b) The fringes have the desired orientation. 

(c) The fringes have the desired spacing. 

(d) The white-light fringes are in the correct location in a 
vertical cross section through the test section. 

(e) The two beams lie in a single plane. 


Evaluation of Density Fields 

The theory of the evaluation of interfere grams of one- and two- 
dimensional flow fields given here is no different from that given elsewhere 
(See reference 11.) The theory is repeated here for the sake of convenience 
The method of measuring fringe shifts is believed to be somewhat different. 

The theory of the formation of the fringes has been discussed in a 
previous section. The production of fringe shifts is illustrated in 
figure. 9* Consider that when no disturbance is present both beams travel 
in air of density p and refractive index n. The value of the refractive 
index is a function of the density of the air, for a given wave length of 
liggit, according to the Lorentz -Lorenz relation 


n 2 - 1 


+ 2 


« p 


or 


n + 1 


(n - 1) 


cc 


p 


n 2 + 2 
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Because n is very nearly equal to unity, the equation can "be written 


n - 1 = kp (l) 


where the constant of proportionality k is the Gladstone -Dale constant* 
Then 


n - 1 


P 


n* 


- 1 
p' 


or 


n’ - n = (n - l)^ - 1^ (2 ) 


Now let one of the two beams pass,. for the distance L, through air of a 
different density p' and index n* , and let values of density and index 
be constant along the length L, but let them be functions of the vertical 
coordinate y. Then the wavefronts will be distorted, as in figure 9. 

The amount of distortion, or retardation X at any point is a function of 
the velocity of li$it at that point. The time for passage of ligit of 
velocity V* through a distance L is L/V* . In that same time, light 
of velocity V will pass* through a distance L + I. Therefore 


L__ _ L + X 
V' V 


But since 


V_ _ X_ 
V' X’ 


L + X 
X 


L_ 

X* 
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or 


X = ~ (X - X.*) 


But, since 


Vo X,o 

n = T ' = T 


and 


n' 


X’ 


therefore 


X = | (n* - n) 

= L ~ (n’ - n) 
K o 


By similar triangles, it follows that 
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Then, "by use of equation (2), 


Y 

h 


L / 

■ XJ (° - 



or 



i (ho 1 N 

b y L n — lj 


+ 1 


Y 

If p the fringe shift in terms of fringe width, is designated hy S(y), 
and the quantity in parentheses is designated by C, then 


P 


CS(y) + 1 


( 3 ) 


In this equation p'/p is the density ratio between some position in the 
disturbance and the undisturbed air* 

The technique of obtaining the values of the fringe shifts for a 
given cross section of the flow field is simplified by the plotting of 
graphs. An interferogram is taken of the undisturbed fringes. (See, for 
example , fig. 10.) Then interfero grams of the flow field are taken. (See 
fig. 11.) Then enlargements to about 7 diameters are made. A position 
along the horizontal axis is chosen at which the density variation is to 
be determined . At that position on the enlargement of the undisturbed 
fringe pattern, a vertical line is drawn, and the positions of the fringes 
are measured with a scale. (A piece of cross-section paper pasted on the 
picture is useful.) Then the fringe positions are plotted as a function 
of the vertical coordinate y. Such a plot is shown in figure 12. (A 
convenient position, such as the lower edge of the nozzle opening, is 
chosen as the zero of the y-coordinate. This position was located on 
the picture in the following manner: Two small pieces of drill rod were 

inserted into the lower nozzle block. From their actual diameter and their 
measured diameter on the picture, the magnification over actual size was 
determined. From the known distances of these rods to the inside edge of 
the nozzle, the location of the edge could be found on the. picture. Two 
small pointers placed iri the camera, the shadows of which can be seen in 
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the lower part of the interf erograms , served to locate, in the pictures 
that did not include the nozzle end, the horizontal line that would extend 
to the nozzle edge.) (The two pieces of drill rod were also used for 
checking the alinement of the light "beam with the nozzle edge. They are 
placed on opposite sides of the center of the nozzle. If the alinement 
is correct, then when, on an interferogram, the proper distance is 
measured from each rod to give the location of the edge of the nozzle 
opening, the same location is obtained for both measurements.) Then, on 
the enlargement of the interferogram of the flow field, at the same position 
horizontally, the positions of the fringes were measured on a vertic al 
scale. Then obtaining the fringe shift as function of y was simply a 
matter of going to the first plot, figure 12, with a value of y and the 
corresponding fringe number, reading there the undisplaced fringe number 
for the same value of y, and subtracting. For example, suppose that 
fringe number 20 lay at y = 17 in the flow interferogram, and fringe 
number 8.4 lay at y = 17 in the no-flow interferogram; then the fringe 
shift at y — 17 is 11.6 fringe widths. The density ratio between the 
place where y = 17 and the undisturbed air is obtained from equation (3). 
This procedure was carried out for each cross section through the flow 
field where it was decided to obtain the density distribution. 

X 1 

The "constant" C = — n _ ^ was determined for each cross section 

through the jet. The value of X Q was 5170 angstroms. The value 

of n - 1 was determined from equation (l). The value of p, the density 
of room air, was obtained from pressure and temperature measurements. The 
best critical-table value of k is O.UL67 cubic foot per slug for 
5170 angstroms. The width of the nozzle opening was 2.999 inches, and 
this value would be the value of L if there were no e nd -effects. A 
corrected value of L was arbitrarily obtained, as follows: The fringe 

shift was plotted against y from the undisturbed region, throu^i the 
boundary layer, into the free- stream region of the jet. The area under 
the curve was obtained and divided by the fringe shift in the free stream. 
This gave a position at which the same fringe shift would occur if it 
took place abruptly, rather than gradually through the boundary layer. 

The difference between this position and the position of the edge of the 
nozzle was subtracted twice from the actual nozzle width to obtain the 
value of L used in the "constant" C. This method is sufficiently 
accurate for cross sections close to the nozzle. 


Effect of Refraction 

When light passes throu^i a medium in which the index of refraction 
varies in a direction that is perpendicular to the initial direction of 
propagation of the light, the direction of propagation is changed, or 
the li^it is refracted. In the mixing region of the jet considered 
herein, the density of the air varies , and consequently also the index 
of refraction varies. The purpose of the present section is to determine 
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whether the refraction, or tending, of the li$it has a significant effect 
on the mixing-region density distribution obtained from interf erograms . 

Let y be the vertical coordinate through the mixing zone and z the 
horizontal coordinate in the direction of the light beam. Because the air 
density in the mixing region increases as y increases, the light will 
be refracted in the manner indicated qualitatively in the following figure: 



If (ft represents the angle of incidence, then Snell's law states that, at 
any point in the medium. 


sin $ = 


JT 

ci 


where Y is the velocity of the light and c 1 is a constant. Accordingly, 
therefore , 


tan <j> 


T 


sin ^ 





V2 

7FF 


dz 


dy _ / c l 2 
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"but 


1 

V 


n 



therefore 



It is now required to solve this equation in order to find the amount by 
which a ray of light is deviated in the mixing region. First, it is 
necessary to express n as a function of y. In a subsequent section, 
the density distribution through the mixing region is obtained. For the 
present purp° se this distribution is approximated by a linear variation 
that fits the actual variation over a large portion of the mixing region. 
The density of the air at the outside of the mixing zone is taben as 
0.0021+ slug per cubic foot and at the inside edge of the mixing region 
as 1.5 times as great, or O.OO36 slug per cubic foot. The effect of 
refraction is greatest at the place where the density gradient is the 
greatest. This occurs at the cross section that is closest to the nozzle. 
For the present investigation that cross section is 2 inches from the 
nozzle. The actual width of the mixing zone there is about 0-33 inch. 

The assumption is made of a linear density gradient equal to the average 
gradient across the mixing region. The assumed density variation, then, 
is given by the equation \| 


p = 0 .002h + 0.0036y 


and the index variation is, by use of equation (l) and the given value 
for k, 


n = 1.0002 8 + 0.00042y 


or 


n = a + by 



The differential equation then becomes 


dz 




By substitution of 


aci 



and 


then 


I = p + qy 


/ 


dz = 


d| 


W~^ 


On integration, 


q(z + 02) = 



For evaluation of the integration constant 
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But at z = y = 0, 


C 1 


= 7 = 




a 


Therefore 


qc 2 = logg 1 = 0 


and 


c 2 = 0 


Therefore 


and 



h 

a 


z = log e 




or 


- z — z / . 

e a + e a =2/l+l 

The path of a ray , therefore, through a medium of linear density gradient 
is a catenary. Because a is very nearly equal to unity, the equation 
can he simplified to 

e^ z + e = 2(l + hy) 
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or, by use Of the series expansion, to 

y £ I bz 2 

The deviation of a ray on passing through the 3 inches of the mixing 
region, from z = 0 to z = 3» is 


y = ~ X 0 .00042 X 9 
= 0.0019 inch 


The index of refraction at y = 0 is 


n = 1.00028 


and at y = 0.0019 is 


n = 1.0002808 


The important quantity, though, is. n — 1 . The light ray emerges from 
the mixing region at a place where n — 1 differs by less than 1 percent 
from its value at the place where the ray entered the mixing region* For 
the jet under discussion, therefore, the effect of refraction is negligible 


RESULTS 
Interfe ro grams 


For obtaining interferograms of the mixing region of the free jet 
of Mach number 1 . 6 , the interferometer was so adjusted that straight, 
horizontal interference fringes were produced when there was no air flow, 
as is shown in figure 10. 

Figure 13 shows the portion of the jet of which interferograms were 
taken for the present investigation. This portion was the bottom part 
of the horizontal, jet for the first 10 inches from the nozzle. The 
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ID -inch length was covered by taking a series of seven interf erograms , each 
of which covered a portion of the jet and its mixing region that was about 

2~ inches high and about inches wide. Only the first six of these 

interferograms were used in obtaining the density and the velocity distri- 
butions that are given in the present paper. These six interf erograms 

covered the first 8^* inches of the mixing region. The interferograms are 

shown as figures 11(a) to 11(f)* A composite of the six interferograms 
is shown as figure 14* 

Three regions are distinguishable in these pictures. The lower portion 
of each picture, where the fringes have retained essentially the same 
spacing that they had with no air flow, was taken with light that passed 
through nearly undisturbed room air. The upper portion of each picture 
shows the free-stream region of the jet. The center portion of each 
picture is the mixing zone between the free room air and the free-stream 
portion of the jet. 


If the irregularities in the fringes in the lower and the upper portions 
of the interferograms are averaged out, then the spacing of the fringes is 
found to remain the same as it was with no air flow. This shows that there 
is uniform density in each of the regions. In the mixing zone, however, 
the fringes are crowded together and there is a density gradient in that 
region. The density varies, in fact, from atmospheric outside the Jet 

to about It- times atmospheric inside the jet. 


As can be seen in figure 11(a), at the end of the nozzle the boundary 
layer that has built up along the nozzle and that emerges from the nozzle 
is much thicker than the beginning of the jet mixing region. At a 
distance of about 2 inches downstream from the nozzle, however, this 
boundary layer has lost its identity and there is only the mixing zone, 
the undisturbed jet, and the room air. Only data takBn from the region 


between 2 and 7“ inches downstream from Idle nozzle are included in the 


present paper. In that region seven vertical cross sections were chosen, 
at 2, 2^, 3^, 4^, 5^, 6, and 7 t; inches from the nozzle. 


Density Distribution 

The density variation along each cross section was obtained. The 
method of obtaining the density variation was first to measure the variation 
of fringe shift along each vertical cross section. This was done by the 
method described in the section entitled "Evaluation of Density Fields ." 
Figure 15 shows a plot of the density variation across the mixing zone at 
the seven cross sections. On the vertical axis is plotted the ratio of 
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density to atmospheric density p/p a tm* the horizontal axis is plotted 
the nondimens ional parameter <ry/x. The variables y and x are position 
coordinates. The y-axis is vertical. The x- axis is not quite 
horizontal hut has been so chosen that it coincides with the line along 
which the density ratio is 1.1. This is' the line along which the velocity 
is 0 .5 free -stream velocity, (it has been customary in the past to place 
the x-axis along the 0.5-velocity-ratio contour. For the jet under 

discussion, this contour was at an angle of -lj^° to the horizontal.) The 

parameter a is an experimentally determined scale factor. Its value is 
obtained by comparing the experimentally determined velocity distribution 
with the theoretical distribution. The parameter o will be discussed 
in a subsequent section. Figure 15 shows that there is fairly good 
similarity in the density distributions. The lack of scatter of the 
experimental points at the outside portion of the mixing region is explained 
by the fact that atmospheric density was vised as the reference in the 
quantity p/ P a tm was plotted. The scatter that occurs at the inside 

part of the mixing region can be attributed partly to variations in 
atmospheric density and partly to variations in the stagnation temperature 
of the Jet , both of which varied from picture to picture . 


Velocity Distribution 

From the density distributions the velocity distributions were 
calculated with -the aid of several assumptions and approximations. It 
was assumed that the static pressure in the jet and the ml -ring region was 
the same as the pressure in the room air outside the jet. The temperature 
distributions through the jet were then obtained from the density distri- 
butions by the general gas law. 

For the different interfero grams the room-air temperature varied 
between 79*6° F and 8o.l° F. The stagnation temperature of the jet air 
varied between 71.6° F and 73*6° F. (The stagnation temperature was 
measured with a thermocouple that was installed in the 6- inch pipe ahead 
of th9 nozzle where the air velocity was low, about 150 feet per second, 
and the temperature recovery factor was very nearly equal to unity.) The 
assumption was made that there was no heat transfer in the mixing region. 
The jet and the mixing region were therefore considered to be isothermal 
from the standpoint of stagnation temperature. For the calculations of 
velocity distribution, the assumption was made that the stagnation 
temperature in the mixing zone was the same as the temperature of the 
room air. Then, from the constant stagnation temperature and the static 
temperature distributions, the velocity distributions were calculated 
from the conservation-of -energy equation 



where u is the velocity and Cp is the specific heat at constant 
pressure. 
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Tollmien obtained (reference l) the theoretical velocity distribution 
for the nri vlng region of an incompressible jet. His results, given in 
table I of reference 1> are shown in figure 1 6. The vertical coordinate 
is u/uq , where Uq is the free-stream velocity. The horizontal 
coor din ate is again ay/x. (For this figure only, the position of the 

x-axis has not been adjusted to coincide with the ~ = 0.5 contour but 

is horizontal, at right angles to the y-axis.) Tollmien's results show 
that the outside edge of the mixing region of an incompressible jet is 

given by = -2.04. Furthermore, Abramovich found (reference 3) that 

compressibility had no effect on the value of c ry/x at the outside boundary 
of the -ring region. He treated theoretically a Jet in which the 
stagnation temperature was the same as that of the ambient air and in 
which the free-stream velocity was large - up to the velocity of sound. 

He also treated the low- speed jet in which the stagnation temperature was 
different from the temperature of the ambient air. He found that compressi- 
bility effects arising from the high velocity or from the temperature 
difference did not affect the value of ay/x at the outside boundary, 
but that the value for both cases was -2 .04. In the present paper, therefore, 
this value has been accepted as correct. 

Tollmien's results also show that the inside edge of the mixing 

region of an incompressible Jet lies along a value = 0 .98. Abramovich 

found, however, that the value of ay/x at the inside boundary should be 
affected by compressibility considerations. His analysis, however, was 
not extended to supersonic velocities, and no theoretical value for ay/x 
at the inside boundary is available for supersonic flow . 


The quantity a was introduced in the theory of free jet mixing as 
a scale factor. (For example, see reference 1.) ,The theory does not give 
the value of a, but its value is determined by fitting the experimentally 
determined velocity distribution to the theoretical velocity distribution. 
In the present case, the inside boundary of the mixing zone is clearly 
shown on the interf erograms , and the angle that it mates with the free- 
stream direction can be measured and shown to be 3°» Because the 
theoretical value for ay/x at this boundary is not known, the value of a 
cannot be determined from the measured rate of spread of the mixing zone 
into the jet. Furthermore, at the outer edge of the mixing zone where 


= -2.04, the density gradient is extremely small. An interferometer 

is not very sensitive to such small gradients, however, and the outer . 
edge of the mixing zone cannot , therefore , be determined accurately from 
the interf erograms . Shadowgrams that covered an extent of 5 feet along 
the jet were therefore taken. On the shadowgrams the outside boundary 
of the mixing zone appeared to lie at approximately 6°. This angle would 
give a value of a of approximately 20 , 

a ~ 

~ -tan 6° 
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It was felt, however, that the shadowgrams might not accurately Indicate 
the true boundary of the mixing zone , and that the most satisfactory 
method of determining a would he to choose a value that would give the 
best fit of the, experimentally determined velocity distribution with 
Tollmien's theoretically determined velocity distribution over the sub- 
sonic portion of the mixing region. First, Tollmien's curve, shown in 
figure 16, was shifted horizontally by 0.39 in order that the 0.5 value 
of u/uq would lie at the zero value of cry/x, as shown by the solid 
curve of figure 17 . Then the velocity distributions that were calculated 
from the measured density distributions, as has been indicated, were 
plotted, with a taken as 15 and with the x-axis taken along the 

— = 0.5 contour. This value of a was chosen because it makes the 
u° 

data agree with To llmi en's curve between about 0.2 and 0.6 on the vertical 
scale. At values of u/uq smaller than approximately 0.2 the data are 
probably not very accurate. At values of u/uq greater than O.625 the 
flow is supersonic and it is to be expected that compressibility effects 
alter the velocity distribution from that of incompressible flow. In 
fact, figure 17 shows that the inside boundary of the mixing region of 

the supersonic jet is at 2Z = 1.12, compared to 1-37 for the incompressible 

x 

Jet. In the supersonic jet, therefore, the mixing region, spreads into the 
free stream at a slower rate than in an incompressible jet- 

The value of a obtained here can be compared with the values 
obtained by other investigators for incompressible jets. Measurements 
in a large jet at Gottingen (see reference l), in which the free-stream 
velocity was about 100 feet per second, gave a value for a of 11.8. 

Liepmann and Laufer state (reference 8) that Cordes found a value of 
11.95 for a. Liepmann and Laufer found that for their jet, which had 
a free-stream velocity of 59 feet per second, the value of a was 12.0. 

These values are to be compared with the results of the present measurements, 
which gLve a value of 15 for a for a jet of Mach number 1.6. Because a 
is a measure of the rate of spread of the mixing region, the mixing region 
of the present supersonic jet spreads less rapidly than that of incompressible 
jets. The rate of spread of the outside boundary is 12 /l5 that of incom- 
pressible jets. The ratio of the rates of spread of the inside boundary 
is even less, as is shown by the fact that, on figure 17, the inside 
boundary lies at a smaller value of ay/x than the theoretical. 

In the theory of the incompressible Jet, it has been customary to 
assume that the value of the "mixing length" is constant across the 
mixing zone and that the 'value of the mixing length l is proportional 
to the distance from the nozzle, 


l = cx 
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The constant of proportionality c is vl/2<JJ. For the subsonic jets, 
with a - 12 , the proportionality constant c has the value of 0 .017 • 

For the supersonic jet, with c = 15, the value of c is 0.012. By the 
use of the hot-wire method, Liepmann and Laufer have shown that the "mixing 
length" is not constant across the mixing zone. If a is taken, notwith- 
standing this fact, to "be a measure of the amount of turbulence, then the 
present measurements show that the mixing region of a supersonic jet is 
less turbulent than the mixing region of a low-speed Jet. 

Figure 17 shows that the velocity distributions at the various cross 
sections are similar. The figure also shows that the flow in the ml ring 
zone is turbulent , because the rate of spread is linear, in agreement 
with theory, as compared to laminar flow, in which the rate of spread is 
a function of the Reynolds number based on the x- coordinate. 

Figures 15 and 17 show density and velocity distributions, respectively 
at seven cross sections. These cross sections lie between 2 and 7i inches 

1 2 

downstream from the nozzle. Interfere grams of the region from 7= to 

10 inches downstream were also taken and density and velocity distributions 
were obtained. These distributions were very similar to those shown in 
figures 15 and 17* The distributions so obtained, however, at these 
distances from the nozzle are not necessarily accurate. At 10 inches, the 
path length of the light through the mixing regions on the two sides of 
the Jet is about 3 inches. The effect of this on the results was reduced 
by the method described in a previous section. When the effect is large, 
however, it is not certain that the method of correcting is at all accurate. 
For that reason no results are shown in the present paper for distances 

greater than 7^ inches from the nozzle. 


SIMMARY OF RESULTS 


It has been found that, for the free supersonic jet of Mach 
number 1.6, 

1. Density distributions through the mixing region were similar to 
each other at the cross sections investigated. 

2 . Telocity distributions through the mixing region were similar 
to each other at the cross sections investigated and were similar to 
Tollmien's theoretical velocity distribution in the subsonic portion 
of the mixing region. 
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3 . The turbulence in the mixing region was less than that for incom- 
pressible jets. 

4. The rates of spread of the mixing region, both into the Jet and into 
the ambient air, were less than those of incompressible Jets. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., January 21, 19 U 9 
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Figure 1 .— Supersonic nozzles 
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Figure 2.- Basic arrangement of Macb-Zehnder interferometer. 
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Figure 4.— Region in which fringes appear to he formed in Machr-Zehnder interferometer. 
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Figure 5.— Full view of interferometer with light source (left) and camera (right). 
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Figure 6 .— Interferometer. (Splitter plates, lower left and upper right 
mirrors, upper left and lower right.) 
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Figure 7 »— Interferometer mirror mount 
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Figure 8 .- Complete optical diagram of interferometer used in this investigation. L x , parabolic 
mirror f/9-0, 36 -inch focal length; L 2> Eastman lens f/ 3 . 5 , 13* 5-inch focal length; 

Eastman lens f/l. 6 , 50-millimeter focal length; Bausch and Lomb f/2.3, 2— inch focal 
length; L^ f Bausch and Lomb f/8.0, 6 — inch focal length; and Co , compensating plates. 
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(a) 0 to 1^ Inches from nozzle. 
Figure 11.— Interferogram. of Jet. 
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Figure 11 .— Continued. 
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(c) 3 to inches from nozzle. 


Figure 11 .— Continued. 
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Figure 11.— Continued. 
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Figure 11.— Continued. 
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(f ) 7 to 8|- Inches from nozzle. 


Figure 11.— Concluded. 
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Figure 12.— Location of undisturbed fringes. 
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Figure IT.- Variation of velocity ratio through mixing region 





